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ABSTRACT 


Gamma ray transport and shielding properties of Lead (Pb) and some multilayers and their reversed combinations have 


°°2Cf point source. Nal (T1) scintillation detector has been used for the detection. The 


been investigated using a 3.7 x 10° Bq 
unfolded dose spectra of the fission gamma ray spectrum for some specific shield thicknesses were also studied. Three 
dimensional Monte Carlo computer code MCNP have been used for theoretical calculations to study the characteristics of 
the above mentioned shielding setups for fission gamma rays. Experimentally obtained values are then compared, 
effectiveness of the materials are satisfactory. It has been found that the materials (single layer and multilayer) studied are 


very effective as the shield for gamma rays. 
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INTRODUCTION 


All penetrating radiations from radioactive substances and their harmful effects has been a matter of concern and study. 
The most significant radiation fields are the combinations of [1], fast neutrons, thermal neutrons, primary gamma rays, and 
secondary gamma rays. In shielding design an important factor that is required to be considered is the choice of most 
suitable material from which the main part of the shield should be built. It also satisfies all other requirements like the 
density, the weight, and the amount of the material needed for shielding purpose. The primary gamma rays is the main 
radiation of interest which can best be shielded with lead or other high density materials [2]. And the last one - the 


secondary gamma rays with an energy of 2.2 MeV are produced as a result of the capture of thermal neutrons by hydrogen. 


The combination of the suitable component shielding materials together is required for increasing the 
effectiveness of the shield. The combination of shielding components into a most effective single homogeneous composite 
has now become a recognized standard method of developing shielding materials. The composites so developed are given 
variety of shapes or forms such as bricks, blocks, slabs, sheets, cylinders, pellets, casks etc. according to the requirements 
[1]. An alternative method is to form suitable multilayers. Such type of heterogeneous mixers [3] can also be developed by 


forming multilayers of the composite shielding materials. 
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In shielding gamma rays traditionally theoretical as well as experimental investigations have been considered. The 
theoretical approach considers the interaction cross section of different materials with a view to developing the required 
shield specifications. The experimental method generally yields gross results and is concerned only with the total effect on 


half value thicknesses and build up factors. 


In reviewing the present trend of studies, it is found that among all the applied radiation sources, diagnostic x-ray 
machines got the priority [4]. Shielding experiments have been performed on a few shields; the results are then compared with 


those obtained through calculations using MCNP code for neutron and gamma ray transmission. The gamma-ray spectrum is 


252 


obtained in the present work from a ~“Cf source. It was necessary to identify the components of the gamma rays, and a part of the 


present work is devoted to this aspect of unfolding study of gamma rays as it simultaneously emits neutrons and gamma rays and 


the source have different energies and intensities. The total dose due to 25, 50 and 100 keV gamma rays unfolded was evaluated. 
THE CONCEPT OF UNFOLDING 
The differential pulse height spectrum N(H) observed from any radiation detector is given by. 

N(H) = JR(H,E)S(E)dE (1) 


Where, R(H,E) is the response function of the detector, and S(E) is the energy distribution of the incident 
radiation. The response function usually is dependent on [5] (i) the operating conditions of the detector, (11) the source 


detector geometry, and (111) the detector counting rate. 
We have a different discrete form of the recorded count in the channel (N;) in the case of a multichannel analyzer, 
N; = VRS; (2) 


Here, Rj; is the response matrix coupling the i pulse height interval with the jth energy interval. S; is the radiation 


intensity in the jth energy interval. 


At times we become interested in deducing the energy distribution of the incident radiation. As the recorded 
spectrum is made up of M values of N; we may write M simultaneous equations like the equation (2). Assuming that the 
response matrix Rj is known, the process solving the equation for all the S; elements for M => L is known as spectrum 
unfolding. Unfolding methods have been thoroughly developed for application area, of spectra recorded by gamma ray 


spectrometers employing either Nal(T1) scientillators or Ge(Li) detectors. 
Unfolding Method 
e Mathematical Formulation 


The conversion is achieved in accordance with the following relation: 


S(E)) AE; = zi * VINE) AE, (3) 
Here N(E;) = The number of counts per sec unit pulse height in bin j 

Ri =The inverted response elements 

J =The maximum energy bin converted 


Now from the equation (2) the operation 
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S=R'N (4) 


Should give the required photon strength, S; but it has a number of disadvantages: (1) direct matrix inversion 
amplifies small statistical fluctuations resulting in violent oscillations and (11) sometimes negative strengths of photons in 


the unfolded spectrum [6] is observed. These disadvantages are avoided by adopting Gold’s iterative technique [7]. 
e The Dose Rate 


The dose rate D now can be expressed as [8] 


poe | 
D = *£i-0 (Ei) EiSCE)) U, (Ei) dE; (5) 


Where, K(E;) = a conversion factor, 

E,; =Photon energy, 

S(E;) = Number of photons calculated per unit energy, per sec and per cm” 
u,(E;) = energy absorption coefficient, defined as 


Ei 
u,(E;) ="! (BE) (6) 


In terms of K; and E; the dose rate may be expressed as 


| 
D =" fi-0 (E;) E;S(E;) w(E\)dE; (7) 


Here E; = average energy transferred by photons per collision in the medium, and 
u(E;) = linear absorption coefficient. 
e Calculation for R 


In the work the response matrix described in Eqn. (1) has been calculated. The following procedure was adopted for the 


purpose. The detector used was a 3"x3" NalI(T1) detector with the following amplifier gains: 


Coarse gain = 4 
Fine gain =. 
Applied high voltage = 750 volts. 


The amplifier used was the model 2012 with positive input polarity. 


The output obtained from the Nal detector was fed to the input of a computerized multi channel analyzer (MCA). 
The MCA was calibrated using some known sources. Eight different known gamma sources (e.g., “41 Am, }°Cd, “Mn, 
'Co, ‘Ba, ”Na, '°’Cs, °°Co) have been used for the purpose. Each of the sources were kept at a distance of 5 to 10 cm 
from the detector. The channels were selected before calibration and times were also preset. After this, the spectrum was 
collected for sources of different energy and source strengths keeping the distance between the source and the detector 5 to 
10 cm and for a preset time. The spectrum was then stored in the hard disc of the computer. The energies together with the 


yield per cent of the sources [9] and their distances from the detector has been shown in Table 1. 
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The strength of the sources on the date of collecting the spectrum was calculated by using radiation decay 


formula. The number of gamma photons incident on the detector surface was found using the following formula: 
S(E) = (A/41)x(a/t’) (8) 


Where, A = Activity of the source on the date of experiment, (a/r*) = solid angle subtended at the source point by 


the detector, a = surface area of the detector = 2d* (d = detector radius), and r = source to detector distance. 
Therefore, S(E) = number of photons incident per second on the detector = (A/4) x (d/r)*. 


If the source is a multiple energy one then the number of photons incident on the detector per second for each 
energy had to be determined. For this the yield (%) had to be added up for all the energies. S(E) was then divided by the 
sum of these yield (%) and subsequently multiplied by the yield (%) of the expected energy. The sum of counts of 10 or 


20 channels (sum count) was divided by S(E) x total time. This gives the efficiency of the detector per photon. 
1.e. Efficiency = [Sum count] / [S(E) x time]. 


The channel number varies from source to source and for variation of their energy as well upto which the 
efficiency has to be determined. The final channel has to be considered as the one upto which the peak area starts. Then 
the efficiency versus energy curves are to be drawn listing different source energies for different range and efficiency. And 
from this, efficiencies are to be calculated for both energies and also for all cases. A matrix has to be made for efficiency 


versus energy from these calculations. 
THE EXPERIMENT 


The experiment was done at Atomic Energy Research Establishment, Savar, Bangladesh. The source used in the work was 
a °°Cf source of strength measured on the date of performance of the experiment. Gamma particles coming out of the 
source was detected by a Nal(T1) detector placed at 105 cm away from the source. Different shielding arrangements were 
placed in between the source and the detector. Spectra were then taken for different shields and different thickness values 


for a preset time. 

RESULTS AND DISCUSSIONS 

Unfolding of Gamma Rays: Unfolded Dose Spectra for the Single Shields 
e PB 


In Fig. 1 all the four different thicknesses for unfolded dose spectra of fission gamma of the shield shows the same 
positions of the peaks between 150 to 450 keV. With increasing shield thickness the peak strength decreases. On the higher 
energy side the dose rate falls at a quicker rate and practically maximum 850 keV gammas also contributes to the total 
dose. On the lower energy side the extension occurs down to 50 keV. The dose rate between 50 to 150 keV is 
approximately one-third of that in between 150 to 250 keV. The total dose values for 4, 16.5, 33.8 and 54 cm shield 
thicknesses are respectively 12.13, 10.39, 7.58, 2.03 uSv/h. Compared to these 100 keV bin results the corresponding total 
dose values found for the 50 keV and 25 keV bins are observed to be the same approximately. The values are given in the 
legends. The peak positions of the main peak of the latter two cases (Figs. 2 and 3) are again found to be the same. This 
means that between 250 to 350 keV maximum numbers of gamma photons are found to possess energy. When the shield 


thicknesses are increased these energy limits does not changes. 
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e Pb 


Figs. 4-6 show unfolded dose spectra of fission gamma after passing through the 5, 10.12, 19.66 and 24.57 cm shield 
thicknesses of Pb. Pb, being a very good gamma shield, reproduces suppressed peak of height 1.2 uSv/h for the 5 cm shield 
thickness for the 100 keV bin spectra. The peak position of the 50 keV bin shifts by about 25 keV towards the lower 


energy side. 


No extra feature (Fig. 6) between 25 to 75 keV in the 50 keV bin spectra could be observed in the 100 keV bin 
spectra. As of the previous case the total dose values are found to remain approximately the same in all the different bin 


spectra. The suppression of the spectra for Pb is due to the high quality shielding effect of the shield. 
Unfolding of Multilayers Containing PB 
e Unfolded Dose Spectra for the Shields PB+Pb and Pb+PB 


The presence of Pb in the multilayers PB+Pb and Pb+PB influences strongly (Figs. 6-7). Hence much depressed spectra 
than those of other multilayers could be observed. The peak position of the main peak is found to be the same as of the 
previous cases. On reversion and thus forming Pb+PB setup less values of dose rate and the total dose rate are obtained 
(please see the legends). The maximum energy gamma photons although very less in number are found to possess 850 keV 


energy. The common points mentioned previously are also found for the PB+Pb and Pb+PB multilayers. 
CONCLUSIONS 


In the present study five single layer shields e.g., PB and Pb and their multilayers e.g., PB+Pb, and Pb+PB and their 
reversed combination have been studied. Flux spectra and dose spectra for fission gamma rays shielding have been 
investigated for all the above setups. The fission gamma ray flux spectra of the PB-multilayer shields shows some special 
features. These features include an initial peak at 0.3 MeV, a trough at 1 MeV, a second small peak observed at 1.5 MeV, a 
trough at 3 MeV, a third peak at 4 MeV and a constant photon flux between 5 to 10 MeV (Figs. 1-3). In all the shielding 
arrangements it is observed that the transmitted gamma photons unfolded by using three energy bins, contain maximum 
number of photons in the range between 250 to 350 keV. On reversing 1.e., by forming a multilayer placing the heavier 
element first, the dose rate values are found to decrease. Beyond the 850 keV energy, practically no gamma photons were 
observed (Figs. 1-14). A multilayer containing more than two layers and having a heavy shield placed first is found to be a 
better gamma shield compared to its reversed setup. Multilayers are found to be more effective at larger thicknesses of the 
shields. Characteristics of PB-multilayers for fission gamma rays are found to be dependent on their uncommon component 
shields. On reversing the components the shielding characteristics change, but the relative characteristics of the uncommon 


components do not. 


Among the single layers studied Pb is the best gamma ray shield. For practical purpose Pb is generally avoided 
because of its toxic effect and then BC naturally becomes the next choice for the whole range of thicknesses studied. For 
fission gamma ray shielding, for small thickness, Pb+PB multilayer shows equal shielding ability as a better choice for 
gamma ray shielding purpose and this is supported by experiments as well as theoretical ( MCNP ) calculations (Table 2). 
Once the choice is made about the suitable shielding arrangement, the results obtained from the present study in the form 
of polynomial fits can give the quantitative information about the exact value of the shield thickness necessary. This is one 


of the beauties of the present investigation. 
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TABLES 


Table 1: Half life, Source Strength, Energy, Yield Percent and The Distance of the Sources from the Detector 
used for the Unfolding =a 


Isotope (T 4/2) srene (keV) (%) Distance (cm) 
eee ee 
Pr osam eertin| [me fs 
jive | stay | gasrsey | 8S ftom | 


years (on 1-7-1983) 356 62.2% 
Oe gig | fae 
Fe foe |i [em 
Table 2: Gamma Half Value, 25% Transmission, Tenth Value and 5% Transmission Fission oe 


Shielding 50% Transmission 25% Transmission 10% Transmission 
Material | MCNP MCNP MCNP MCNP 
[pp | 44 s30_ | ao | Ps PP 





| Pb | 25 | BT OS O83 8S 
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Figure 1: Unfolded dose Spectra (100 keV bin) of Fission Gamma after Passing through 
different Thicknesses of the Shield PB. 
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Figure 2: Unfolded dose spectra (50 keV bin) of Fission Gamma after Passing through 
different Thicknesses of the Shield PB. 
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Figure 3: Unfolded Dose Spectra (25keV bin) of Fission Gamma after Passing through different 
thicknesses of the shield PB. 
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Figure 4: Unfolded Dose Spectra (100keV bin) of Fission Gamma after Passing through different 
Thicknesses of the Shield Pb. 
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Figure 5: Unfolded Dose Spectra (50keV bin) of Fission Gamma after Passing through different 
Thicknesses of the Shield Pb. 
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Figure 6: Unfolded Dose Spectra (25keV bin) of Fission Gamma after passing through different 
Thicknesses of the Shield Pb. 
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Figure 7: Unfolded Dose Spectra (100keV bin) of Fission Gamma after Passing through different 
Thicknesses of the Shield PB+HC. 
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Figure 8: Unfolded Dose Spectra (50keV bin) of Fission Gamma after Passing through different 
Thicknesses of the Shield PB+HC. 
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Figure 9: Unfolded Dose Spectra (100keV bin) of Fission Gamma after Passing through different 
Thicknesses of the Shield PB+Pb. 


Wwww.laset.us editor@tiaset.us 


102 M. S. Hossain, S. M. A. Islam, M. N. Islam, M. B. Ullah & Meherun Nahar 





eT —O— PB=dem+Pb=2.8em,TD=3,121 wSvh 
2 12 } i = | | 
* —O— PBe12.2cm*Pb=2.8em, TD= 1.30018 psv/h 
= 10 
3 —A— PB 12.2em+Pb=8.3em,TD=.4341 wSvih 
= Hi 
= «¢] —¥— PB=24,9em+ Pb=8.3em,TD=0.02245 pSv/h 
rt) , o— 
4 
ao (4 —K— PB=24.9enr+ Ph 18. Gem, TD 0,004245 Swh 
‘ a a _ 
0 


0 100 200 200 400 900 600 foo 600 900 #1000 1100 1200 


Energy in keV 


Figure 10: Unfolded Dose Spectra (S0keV bin) of Fission Gamma after Passing through different 
Thicknesses of the Shield PB+Pb. 
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Figure 11: Unfolded Dose Spectra (25keV bin) of Fission Gamma after Passing through different 
Thicknesses of the Shield PB+Pb. 
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Figure 12: Unfolded Dose Spectra (100keV bin) of Fission Gamma after Passing through different 
Thicknesses of the Shield Pb+PB. 
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Figure 13: Unfolded Dose Spectra (S0keV bin) of Fission Gamma after Passing through different 
Thicknesses of the Shield Pb+PB. 
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Figure 14: Unfolded Dose Spectra (25keV bin) of Fission Gamma after Passing through different 
Thicknesses of the Shield Pb+PB. 
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